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In supersymmetric theories where the lightest supersymmetric particle is the gravitino the next to 
lightest supersymmetric particle is typically a long lived charged slepton. In this paper, following our 
earlier proposal [H, we perform a detailed study of the production of pairs of these particles induced 
by the interactions of high energy cosmic neutrinos with nucleons in the earth, their propagation 
through the earth and finally their detection in neutrino telescopes. We investigate the charged 
slepton energy loss in detail and establish that the relatively small cross-section for the production of 
supersymmetric particles is partially compensated for by the very long range of these heavy particles. 
The signal, consisting of two parallel charged tracks emerging from the earth, is characterized by a 
track separation of a few hundred meters. We perform a careful analysis of the main background, 
coming from direct di-muon production, and show that it can be separated from the signal due to its 
characteristically smaller track separation. We conclude that neutrino telescopes will complement 
collider searches in the determination of the supersymmetry breaking scale, and may even provide 
the first evidence for supersymmetry at the weak scale. 



PACS numbers: 11.30.pb, 13.15+g, 12.60.jv, 95.30.Cq 
I. INTRODUCTION 

One of the most pressing questions in particle physics 
is the origin and stability of the hierarchy between the 
weak and the Planck energy scales. Natural solutions of 
this so called hierarchy problem require new physics at 
the TeV scale. One of the most attractive candidate the- 
ories is weak scale supersymmetry. Although this is in no 
small measure due to its theoretical appeal (it is a simple 
and natural extension of the usual space-time symme- 
tries), it is also favored by data from electroweak observ- 
ables. These point to a weakly coupled Higgs sector, one 
without significant deviations from the standard model 
in regard to electroweak precision observables, as arises 
naturally in supersymmetric theories. However, super- 
symmetry must be broken since the superpartners have 
not yet been observed. The supersymmetric spectrum is 
determined by the supersymmetry breaking mechanism. 

Supersymmetric theories typically possess a discrete 
symmetry, R-parity, which ensures that corrections to 
precision electroweak observables are small. The exis- 
tence of this discrete symmetry immediately implies that 
the Lightest Supersymmetric Particle (LSP) is stable. 
Which of the supersymmetric particles is the LSP is de- 
termined by the scale of supersymmetry breaking, which 
we denote by y/F. When supersymmetry is broken at 
high scales such that \/F > 10 10 GeV the LSP is typically 
the neutralino. If however supersymmetry is broken at 
lower scales, \/F < 10 10 GeV, the LSP tends to be the 
gravitino. In models where the LSP is the gravitino, the 
Next to Lightest Supersymmetric Particle (NLSP) is usu- 
ally a charged slepton, typically the right-handed stau. 
When the supersymmetry breaking scale \/~F is much 



larger than a TeV the stau NLSP decays to gravitinos 
through interactions that are extremely small, and there- 
fore its lifetime can be very large. In gauge-mediated 
SUSY breaking, for instance, we have 
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where mf R is the stau mass. Thus, for \/F > 10 7 GeV 
if these NLSPs were to be produced by very high energy 
collisions they could travel very long distances before de- 
caying. 

Many interesting and realistic supersymmetric mod- 
els have been proposed where the scale of supersymme- 
try breaking \JF is larger than 5 x 10 6 GeV and where 
the LSP is the gravitino while the NLSP is a long-lived 
stau. These include models of low energy supersymme- 
try breaking such as gauge mediation 0, but also super- 
WIMP scenarios where the supersymmetry breaking 
scale is much higher and the NLSPs decay to the grav- 
itino only on very large time scales of about a year. 

In a recent letter [l| we proposed that the diffuse 
flux of high energy neutrinos colliding with the earth 
could produce pairs of slepton NLSPs which due to their 
large range could travel large distances through the earth 
and be detected in neutrino telescopes. The neutrino- 
nucleon reaction can produce a pair of supersymmetric 
particles that will promptly decay to a pair of NLSPs 
and standard model (SM) particles. This process re- 
sults in NLSPs which typically have a very high boost 
and therefore will not decay inside the earth provided 
the supersymmetry breaking scale y/F > 10 7 GeV. (For 
5 x 10 6 < \/F < 10 7 GeV a significant fraction of the 
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decays will occur inside the earth.) Since the NLSP is 
charged, its upward going tracks could in principle be de- 
tected in large ice or water Cerenkov detectors, such as 
IceCube [|| . This is in analogy with the standard model 
charged current interaction giving muons, the primary 
signal in neutrino telescopes. The high boost and large 
range of the NLSPs implies that the tracks are parallel 
and well separated, enabling them to be distinguished 
from background events. Various aspects of this scenario 
were further considered in Refs.[j| H, 0]- 

In principle one could think that cosmic rays could be 
an even better source of supersymmetric events, and ul- 
timately of NLSPs. After all, if one considers a proton 
primary, squark pair production should be larger than 
the slepton production via neutrino interactions, due to 
the fact that it goes through a larger coupling. How- 
ever, this cannot compete with the soft, forward cross 
section of cosmic rays with air, which is large enough to 
stop primaries in the atmosphere. Typical cross sections 
for cosmic ray primaries in the earth's atmosphere are 
around 10 — 100 mb, so that the atmosphere is several 
interactions lengths. On the other hand, supersymmetric 
production cross sections at the energies of interest are 
in the order of 10 — 100 pb, or 10 -9 times smaller. Thus, 
the probability of producing supersymmetric events with 
proton or other cosmic ray primaries is very small. Neu- 
trinos, on the other hand, will always go through the 
atmosphere, with enough of them interacting inside the 
earth within the range of the NLSP from the detector. 
Then, since the neutrino flux is not expected to be that 
much smaller than the cosmic ray flux, using neutrinos 
is advantageous. 

In this paper we present a detailed analysis of the orig- 
inal proposal [l|. We perform a careful study of the 
production of supersymmetric particles, their propaga- 
tion through the earth, and their detection in neutrino 
telescopes. We compute the expected number of signal 
events in IceCube, their characteristic energy distribu- 
tion and the energy deposited in the detector. We also 
investigate the main background, which arises from di- 
muon production, and verify that it can separated from 
the signal by making suitable cuts. 

In analyzing the signal, it is particularly important to 
carefully incorporate the energy losses involved in the 
propagation of sleptons through the earth. Although 
Bremsstrahlung and pair-production energy loss become 
less important for heavy charged particles, photo- nuclear 
energy loss potentially presents a problem. In fact, at 
first this appears to be a great obstacle since it is known 
that photo-nuclear energy loss dominates the propaga- 
tion of tau leptons [1, Q. However, this is not the case. 
As we already sketch in Ref. [l[ and was shown by the ex- 
plicit calculation in Ref. there is a mass suppression 
of the photo-nuclear energy loss, which manifests itself 
for masses well above that of the r and affects the en- 
ergy loss of a particle such as the Ir. We will reproduce 
the results of Ref. @ regarding the photo-nuclear energy 
loss of a slepton, and explicitly show how the mass sup- 



pression comes about and why it is not operative for the 
t energy loss. As a result the photo-nuclear energy losses 
for heavy charged particles, whether scalars such as the 
NLSP or fermions such as heavy leptons, are suppressed 
due to their large masses. 

The main background to the signal arises from di-muon 
events which also give rise to pairs of tracks in the detec- 
tor. However, we show that these can be readily distin- 
guished from the signal by the separation between tracks. 
The key point is that the NLSP range is much larger than 
the muon range, which is only of order a few kilometers. 
Therefore while the typical track separation in the detec- 
tor for the NLSP signal is of order a few hundred meters, 
direct di-muon production results in tracks that are typ- 
ically closer together by an order of magnitude. 

The plan for the paper is as follows: we review in detail 
the relevant cross section for supersymmetric production 
in Section [TTJ In Section IIHI we study in detail the en- 
ergy loss of sleptons as they travel through the earth. 
An analysis of the signal and backgrounds is presented 
in Section IIVI and we conclude with a discussion of our 
results in Section IVl 



II. CROSS SECTIONS 

In this section we summarize the calculation of the pro- 
duction cross section for the NLSP pair. The supersym- 
metric process of interest involves the t-channel produc- 
tion of a slepton and a squark via gaugino exchange. For 
simplicity, we neglect mixing with Higgsinos in the gaug- 
ino sector. Then the processes analogous to the charged 
current (CC) in the SM is uN — > i^q, where q can be 
an up or down-type squark, and the particle exchanged 
in the t-channel is a chargino. These are shown in Fig- 
ure HJa) and (b), and are the dominant contributions. 
The neutrino, always produced left-handed by the weak 
interactions, can interact either with a left-handed down- 
type quark (a), or with a right-handed up- type anti-quark 
(b). This results in the partonic cross sections: 

da^_ = Ml 
dt 2 sin 4 #w s(t-M|)2 1 ■ ' 

~dT ~ 2 sin 4 <V s 2 (i-M?) 2 ' (,) 

where s, t and u are the usual Mandelstam variables, 
and Mu), and rriq are the chargino, the left-handed 
slepton and the squark masses respectively. The left- 
handed slepton and the squark decay promptly to the 
lighter "right-handed" slepton plus non-supersymmetric 
particles. We also include the subdominant neutralino 
exchange, Figure Q] (c)-(d). 

We take = 250 GeV, = 250 GeV and three 
values for the squark masses: m q — 300, 600 and 
900 GeV. These are very representative values in the sce- 
narios under consideration. This is clearly the case for 



3 



11 



(a) 



w 3 



(c) 



u a 

(b) 

V V 

1 w 3 



(d) 



FIG. 1: Feynman diagrams for supersymmetric particle pro- 
duction in vN collisions Charged current (chargino) interac- 
tions: (a) Left-left interaction requiring the insertion of the 
gaugino mass in the t-channel line, (b) Left-right interaction. 
Neutral current: (c), (d). There are analogous diagrams for 
anti-neutrinos as well as for strange and charm initial quarks. 
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gauge mediated supersymmetry breaking, where super- 
partner masses are generated by their interactions. Then, 
typically the tr is the NLSP, being heavier only than the 
ultra-light and very weakly coupled gravitino. Charginos 
and neutralinos tend to be heavier since they also feel 
the SU (2)l interactions. Finally, squarks are heavier still 
since their masses are affected by the strong interactions. 
In Figure [5] we plot the cross sections for supersymmet- 
ric production in i>N interactions as a function of the 
neutrino energy. Also plotted for comparison is the SM 
charged current (top curve) the di-muon (second from 
top) cross sections. 

Every supersymmetric event produces a 1^ and q as 
primaries, which then decay. The resulting cascades al- 
ways leave a pair of £r and SM particles. In what fol- 
lows we will take mj = 150 GeV. The SUSY cross sec- 
tions are considerably suppressed with respect to the SM, 
even when well above threshold. This is to be expected 
since the SM cross section is dominated by very small 
values of x, whereas for the SUSY processes one always 
needs to be above a rather large threshold, resulting in 
x > mg lJSY /(2MpE^), where msusy is the typical mass 
of the supersymmetric particles being produced. How- 
ever, we will see that this will be compensated by the very 
long range of sleptons compared with the muon range, 
which is only a few kilometers. 



FIG. 2: vN cross sections vs. the energy of the incident neu- 
trino. The three lower curves correspond to m; = 250 GeV, 
m„ = 250GeV; and for squark masses rriq — 300 GeV 
(dashed) , 600 GeV (dot-dashed) and 900 GeV (dotted). The 
top curve corresponds to the SM charged current interactions 
and the middle one to the di-muon background. 



III. NLSP ENERGY LOSS 



Once produced by the vN interactions in the earth, 
the pair of NLSPs should range into the detector, just as 
the muons produced by CC events [l(| ■ Charged particles 
lose energy due to ionization processes as well as through 
radiation. The average energy loss can be written as [IJJ 



- — = a(E) + b(E) E , (3.4) 



where a{E) characterizes the ionization losses, and b(E) 
includes the contributions to radiation energy losses from 
various sources including bremsstrahlung, pair produc- 
tion and photo-nuclear interactions. 
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A. Ionization 

Ionization energy loss can be approximated by the 
Bethe-Block formula [ill 



dE 
dx 



P 2 A a ' 
4.m 2 e c 2 v 2 ^j 2 (fa) 2 



2nN A r 2 e m e c 2 
In 



P 



2f3 2 



-5- 



C_ 

z„ 



(3.5) 



In eqn. (|3.5[) , Na is Avogadro's number, r e is the classi- 
cal electron radius, m e its mass, Z a and A a the atomic 
number and mass of the absorber respectively, and z the 
charge of the incident particle in units of e, the electron 
charge. The mean ionization potential / is typically of 
about 100 eV (136 eV for Standard Rock). The density 
correction 5 reflects the effect of the polarization of the 
medium as the incident particle's energy increases. This 
reduces the energy loss at high energies. For the case of 
the NLSP in question, this reduction is still small and 
of the order of 10%. Finally, the shell correction C tak- 
ing into account the effects of atomic bounding, are only 
important for small values of fa, similar to the bound 
electrons'. For fa ~ 0.3 this correction is about 1%. For 
the much larger values of fa relevant here, this correc- 
tion is negligible. Therefore, the ionization loss can be 
parameterized by 



a(fa) ~ 0.08 



MeV cm 2 
gr 



(17 + 2 In fa), 



(3.6) 



and is rather independent of the particle mass. For in- 
stance, for a density of p — 2.5 gr/cm 3 , roughly corre- 
sponding to standard rock in the earth's crust, and for 
fa = 10 5 , we have a ~ 8 MeV/cm. However, at these 
values of fa most of the energy loss will occur through 
radiation. 



B. Radiation Energy Loss 

At sufficiently large energies the energy loss due to 
radiation processes dominates over ionization for all 
charged particles. For instance, for muons this domi- 
nance takes place at energies of several hundred GeV. 
The radiation loss can be written as 

b(E) = ^J^ y d -f y dy, (3.7) 

where Na is Avogadro's number, A is the atomic mass 
of the target, and y is the fractional energy loss of the 
(s)lepton 



E-E' 
E 



(3.8) 



with E' the final lepton energy in the target's rest frame. 
Although there is no simple scaling with mass, heavier 
particles such as NLSPs will have suppressed radiative 
energy losses. In what follows we discuss the differen- 
tial cross section da/dy for radiative energy loss of NL- 
SPs due to Bremsstrahlung, Pair Production and Photo- 
nuclear interactions. 



1. Bremsstrahlung 

The differential cross section for Bremsstrahlung is re- 
viewed in detail in Reference. [l2j . The contribution from 
the interaction with the screened nucleus dominates and 
is given by 



da\ 

^1/ J brem.,nucl. 



a ( 2Z—r e 



:V + y' 



where mp is the mass of the interacting (s) lepton, 



y 
(3.9) 
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and we have defined 



6 = 



mey 



2E(1 - y) ' 



and 



D n = 1.54 A 



(3.11) 



Finally, r e is the classical electron radius and in this case 
B = 182.7. There is also a subdominant contribution 
from the Bremsstrahlung induced by atomic electrons. 
This is down by a factor of Z with respect to eqn. (|3.9[) . 
It is possible to approximate this effect by replacing Z 2 
in eqn. (|3.9p by Z(Z + 1). In order to obtain 6b rem .(S) 
we must integrate y, with y m i n = and y max = 1 — 
(3/4)yeZV3 (m t /E). 

As it can be seen by comparing Fig. Q and Fig. (f5]), 
the Bremsstrahlung contributions are greatly suppressed 
for the slepton case compared to the muon. This can 
already be seen in the expression (|3.9[) . by the appearance 



of a factor 1/m 2 . Thus, for the NLSP case, with masses in 
the hundreds of GeV, this contribution will be negligible. 



2. Pair Production 

As was noted in Ref. @ , the dominant contributions to 
energy loss due to pair production come from diagrams 
which are independent of the spin of the heavy lepton. 
Then it is possible to accurately estimate this contribu- 
tion to b(E) by considering a heavy lepton, instead of a 



5 



0.2 



0.05 



: ■- 

- ' 


jair 












: - 


jrem 
ph nu 


c 










: 

























































































io 4 io 5 io 6 io 7 io 8 io 9 



(GeV) 



FIG. 3: The muon energy loss function b(E) defined by 
eqns. (f3~4)) and ((3~7)l . 
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FIG. 4: The scattering l(k)N{P) 
elastic picture. 



£(k')X in the deep in- 



150 GeV in Fig. ([5|) (dot-dashed curves). Once again, it is 
clear that the energy loss is suppressed in the case of the 
heavier (s)lepton. Although pair creation is an important 
source of energy loss in the slepton propagation, it is not 
the most important one, as we will see next. 



scalar lepton as it would be the case for the NLSP. The 
fermion and scalar cases are almost identical at the high 
energies considered here. A compete QED calculation of 
the direct pair production was done in Ref. [13| . Based on 
this, a useful parametrization for the double differential 
cross section is given in Ref. [lij : 



d 2 a 
dydp 



i-y 

v 



(3.12) 



where A e is the electron's Compton wavelength and p = 
(E + — E~)/(E + + E~) is the asymmetry parameter of 
the pair. The functions </> e and </> M correspond to different 
QED diagrams and include corrections for atomic and nu- 
clear form-factors. Their form can be found, for instance, 
in Reference [l2j]. Just as in the Bremsstrahlung case, 
and in order to take into account the contribution from 
atomic electrons, we can simply replace Z 2 in eqn. (|3.12[) 
by Z(Z+l). Then, the energy loss due to pair production 
is computed as 



h (i?\- Na 



, \ y-r—r-dydp. (3.13) 

A Jy min Jo d ydp 



Here, the integration limits are 
4m 



E ~ y ~ 4 V E 
and 

0< p <(l 



(3.14) 



(l-y)E 2 



1 - 



4m e 
~y~E 



The energy loss by e + — e pair production for a muon 
is plotted in Fig. ([3]) and that for a slepton with m» = 



3. Photo-nuclear Energy Loss 

For a heavy particle such as the slepton NLSP, the 
most important source of energy loss is the photo-nuclear 
interactions. These already dominate the r energy loss. 
However, as we will see below, there will still be a con- 
siderable suppression due to the NLSP mass. This mass 
dependence is subtle, so we spend some time discussing 
its origin. We compute the scattering of a scalar such 
a the t NLSP off a nucleon following the deep inelastic 
formalism. In general the cross section can be expressed 
as 



da 



1 



a 



EE'Q' 



(3.15) 



where E(E') is the incoming (outgoing) electron energy, 
k' is the outgoing slepton momentum, and Q 2 = —q 2 = 
(k — k 1 ) 2 is the photon virtuality, k and k' are the four- 
momenta of the incoming and outgoing slepton respec- 
tively. The momentum fraction of the struck quark is 
described by the Bjorken variable 



Q 2 



Q^ = 

2P.q 2M(E - E') 



(3.16) 



with P and M the four-momentum and mass of the nu- 
cleon respectively. The hadronic tensor is parametrized 
as usual in terms of two structure functions, F\(x,Q 2 ) 
&nc\F 2 (x,Q 2 ): 



wv = 



- M 9u„ 



F2 P P 
M 2 Ey 



F 2 P^q„ + Pu% 
M q 2 



M 2 q 2 ) q 2 
(3.17) 



6 



In the Callan-Gross limit, the hadronic form-factors are 
related by _ 10 - 7 



2xF 1 {x,Q 2 ) 



AM 2 x 2 s 



More generally, departure from this limit, can be 
parametrized by the function R(x,Q 2 ) defined by 



R{x, Q 2 



Fl(x,Q 2 ) 



2xF 1 (x,Q 2 ) ' 
where the longitudinal structure function is 

AM 2 x 2 s 



(3.19) 



F L {x,Q 2 



1 



Q 2 



F 2 {x 1 Q 2 )-2xF 1 (x,Q 2 



(3.20) 

Thus, we can eliminate Fi(x,Q 2 ) in favor of F2(x,Q 2 ) 
and R(x,Q 2 ), where the latter function vanishes in the 
deep inelastic limit. 

On the other, hand the slepton tensor is given by 



= (fc + *')/.(* + *')* 



(3.21) 



We can now compute the double differential cross section 
in terms of the structure function F 2 (x, Q 2 ) as well as the 
Callan-Gross violating function R(x,Q 2 ). The result is 
given by 

d 2 a _ na 2 F 2 (x,Q 2 ) ^ 2 ^ 2 



dydQ 



v 



Y i 



Ami 



1 



AM* 



(l + R(x,Q*)) 



(3.22) 



which agrees with Ref. @ • In order to obtain the photo- 
nuclear energy loss function we have 



h (T?\ NA 

Oph-nuc(,-C/J — 



dy 



da 



dydQ 2 



(3.23) 



{(M+m w ) 2 -M 2 )/2ME, y max = 1-m-jE, 
M 2 ) and the minimum 



where y m 

Qma X = ZMEy - ((M + m,f 
photon virtuality is given by 



Q 2 ain = 2[E 2 (l-y)-kk' ~m 2 ] 



2 2 

(i-y) ' 



(3.24) 



and corresponds to a forward slepton. The last expres- 
sion in eqn. (l3.24p reflects the approximation for small val- 
ues of the slepton energy loss (y <C 1). The Callan-Gross 
limit is an excellent approximation in the energies of in- 
terest here, so we will take R(x, Q 2 ) = from now on. 
We must now specify the structure function F 2 (x,Q 2 ). 
We use the ALLM parametrization [l5], [lj|, which fits 
nicely all data from very low Q 2 to Q 2 = 5000 GeV 2 . 
The parametrization has a x 2 /ndf = 0.97 for a total of 
1356 points. The details of the parametrization used here 
can be found in the update of Ref. [l(| • 
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FIG. 5: The slepton energy loss function b(E) defined by 
eqns. pi} and 



In addition, we must account for nuclear shadowing, 
which can be an important effect at low values of x. This 
is done by defining the ratio 



a(A,x,Q 2 ) = 



F 2 A (x,Q 2 ) 



AF^i 



(3.25) 



where F 2 A (x,Q 2 )/A and F 2 N (x, Q 



2 are the nuclear and 
the free nucleon structure functions respectively. A 
parametrization of this ratio describing the data is given 
in Ref. |8|. It is found to be rather Q 2 -independent and 
is given by 



A -o.i 

„[A..c) = { J 4 - 0691 °sio+ - 097 



x < 0.0014 
0.0014 < x < 0.04 
0.04 < x 



Assuming that Z = A/2, the structure function F 2 can 
be obtained as 

F 2 A =a(A,x)j(l + P(x))F 2 ° . (3.26) 

In ean. (|3.26p F 2 is the proton structure function and 
P(x) = 1 - 1.85x + 2.45a; 2 - 2.35x 3 + x 4 describes the 
Fg/Fg ratio as a fit to the BCDMS data. 

The photo-nuclear energy loss can now be computed 
using eqn. (|3.23|) with the structure function F 2 (x,Q 2 ) 
that is extracted from data and parametrized by ALLM. 
We plot the result in Fig. [5] (solid curve). 

We observe that the photo-nuclear energy loss appears 
to be also suppressed for higher masses when compared 
with the muon case, as can be seen from Fig. [3] (solid 
curve). This seems odd at first, since the effect is larger 
for the r than it is for the /i Q . In order to under- 
stand this mass dependence for the heavier masses, we 
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should note that the cross section that determines the 
structure function F%(x, Q 2 ) is dominated by physics at 
rather low Q 2 (~ 1 GeV 2 ), where resonances and other 
non-perturbative effects set in. However, we see from 
eqn. (|3.24p that for large slepton masses the minimum 
value of the photon virtuality is typically larger than 
this, unless the energy loss y is very small. Thus, the 
most important contributions to the structure function 
are left out as a consequence, leading to a suppression 
for sufficiently large masses m| S> lGeV 2 . This explains 
why the effect is still large for the r but then is suppressed 
for heavier (s)leptons. 

Finally, we consider the possible electroweak contri- 
butions to the slepton energy loss, where the photon is 
replaced by either a gauge boson or a gaugino. Since 
the propagating slepton is "right-handed", it only has 
neutral interactions. Of these, the dominant contribu- 
tion comes from the exchange of a Z boson. For this 
to begin to be comparable to the photo-nuclear process 
described above, we need Q 2 <, M|. However, as men- 
tioned earlier, the structure functions are dominated by 
the low Q 2 region. Thus, we expect the Z-exchange con- 
tribution to be suppressed by this effect. The results of 
the explicit calculation confirms this expectation. The 
energy loss due to Z-exchange is at least two orders of 
magnitude smaller than the photo-nuclear energy loss, 
and even smaller for most of the energies of interest. As 
a source of energy loss it is also smaller than the pair- 
production, and is only larger than bremsstrahlung for 
slepton energies above 10 7 GeV. 

In what follows we will make use of the calculations 
of the energy loss for an NLSP in order to compute the 
number of events ranging into the IceCube detector. 



IV. SIGNALS IN NEUTRINO TELESCOPES 
A. Neutrino Flux 

In order to compute the event rates in neutrino tele- 
scopes, we need to know the incoming neutrino flux. The 
presence of cosmic neutrinos is expected on the basis of 
the existence of high energy cosmic rays. Several esti- 
mates of the neutrino flux are available in the literature. 
In most cases, it is expected that km 3 neutrino telescopes 
will measure this flux. Here we will make use of two cal- 
culations of the neutrino flux in order to compute the 
number of expected NLSP events at neutrino telescopes. 

First, we consider the work of Waxman and Bahcall 
(WB) 17], who pointed out that the observed cosmic 
ray flux implies an upper bound on the high energy as- 
trophysical neutrino flux. This argument requires that 
the sources be optically "thin", meaning that most of 
the protons escape and only a fraction of them interact 
inside the source. In order to determine the neutrino 
spectrum, WB fix the cosmic ray spectral index to —2, 



giving 

(d4u\ _ (l-4)xl0-8 a , 

U^Jwb" & GeVcm S Sr 1; 

(4.27) 

where the range in the coefficient depends on the cosmo- 
logical evolution of the sources. Using the upper end we 
obtain the so called WB limit. 

On the other hand, Manheim, Proterhoe and Rachen 
(MPR) [l8[ obtain and upper limit on the diffuse neutrino 
sources in a way that is somewhat different than WB. In- 
stead of assuming a fixed cosmic ray index for the spec- 
trum, MPR determine the spectrum directly from data 
at each energy. In the case when the sources are con- 
sidered optically "thin" this procedure results in a limit 
similar to that of WB for energies between 10 7 and 10 9 
GeV and larger otherwise. 

We consider an initial flux containing both and 
v e (in a 2 : 1 ratio). Since the initial interactions (see 
Figure [1} produce £l and these are nearly degenerate in 
flavor, the flavor of the initial neutrino does not affect 
our results. For the same reason, the possibility of large 
mixing in the neutrino flux is also innocuous here. 



B. NLSP Signals 

Having the neutrino flux from the previous section, 
the production cross section for NLSP pairs (Section HI)) . 
and their energy loss (Section lllip through earth, we can 
now compute the number of NLSP events at neutrino 
telescopes. In order to correctly take into account the 
propagation of neutrinos and the NLSP in through the 
earth, we make use of a model of the earth density profile 
as detailed in Reference [l9j . 

In Figure [J5] we show the energy distribution for the 
NLSP pair events for three choices of squark masses: 
300 GeV, 600 GeV and 900 GeV. Also shown are the 
neutrino flux at earth in the WB limit, as well as the en- 
ergy distribution of upgoing /i's, and that of direct 
production (see Section [TV C[) . We see that, even for the 
heavier squarks, it is possible to obtain observable event 
rates. In Table U we show the event rates for pair 
production per year and per km 2 . The rates are given 
for the WB flux as well as for the Mannheim-Protheroe- 
Rachen (MPR) flux [l8|, both for optically thin sources. 
For comparison, we also show the rates for the 
background. Thus, km 3 Cerenkov detectors such as Ice- 
Cube, appear to be sensitive to most of the parameter 
space of interest in scenarios with a relatively long lived 
NLSP. The rates are comparable to the ones we obtained 
in Ref.0. 

The NLSPs arc produced in pairs very far from the de- 
tector and with a very large boost. Although the angle 
between the two NLSP tracks, is small due to the boost, 
the large range means that there will be a significant sep- 
aration between the tracks. The question is then whether 
this separation is large enough to be observed, but not 




E v (GeV) 



FIG. 6: Energy distribution of £r pair events per km 2 , per 
year, at the detector. Curves that do not reach y axis; from 
top to bottom: rriq — 300, 600 and 900 GeV. Here, m- t — 
150 GeV and = 250 GeV. Also shown are the neutrino flux 
at earth and the fi and the di-muon flux through the detector 
(curves that reach y axis; from top to bottom respectively). 
In all cases we make use of the WB limit for the neutrino flux. 



too large so as to be larger than the dimensions of the 
detector. A rough estimate can be made by defining the 
separation between NLSP tracks as 8R ~ L 9, with L the 
distance to the production point, typically a few 1000 Km 
and 6 ~ PgUSY/Pboost — ( 10 ~ 4 _ !0~ 5 ), the typical angle 
between the NLSPs in the lab frame. Then we should ex- 
pect that the typical separation between the two NLSP 
tracks is of the order of SR ~ few x (10 — 100)m. 

We performed a more detailed calculation in order to 
have a more precise prediction for the distribution of the 
track separation of the events. We consider the angular 
distribution of the two primaries in the center of mass, 
as implied by eqns. (|2.2[) and (|2.3|) as well as those of 
the subdominant neutral exchange. We assume that the 
angular distribution of the right handed slepton pair in 
the center of mass is the same of that of the primary 
decay particles, the £l and q. This is a good approxima- 
tion for energies well above threshold, where most events 
come from, as it can be seen in Fig. [5] In this limit, 
it is possible to define a "thrust axis", around which 
the decay products of the two primary final states are 
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TABLE I: Number of events per km 2 per year for different 
neutrino fluxes at the Earth. The la mass is 150 GeV and 
squark masses are 300, 600 and 900 GeV respectively. The sig- 
nal events are integrated from threshold (E v > 1.6xl0 5 GeV), 
whereas the number of di-muon events are given for neutrino 
energies above 10 3 GeV. The column corresponds to the 

di-muon background before track separation cuts are applied 
(see Section IFVC)) . 
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boosted. We then boost this distribution to the labora- 
tory frame. The approximation made fails for energies 
just above threshold, for which decay products can have 
a more spherical distribution resulting, in principle, in 
larger angles and ultimately in larger track separations. 
However, the events among these that do make it to the 
detector tend to come from smaller distances than the 
range, therefore compensating -at least in part- the po- 
tentially larger angles. In any case, this approximation 
tends to underestimate the stau track separation. In the 
calculation we also include the additional zenith angular 
dependence induced by the earth profile. 

Taking the above features into consideration, we per- 
formed a simulation where we generated approximately 
30k events (the numbers for each squark mass varied ac- 
cordingly to the production energy threshold). These 
events where distributed in energy steps ranging from the 
stau production threshold to 10 11 eV. For each energy we 
generated events for which the center of mass (CM) an- 
gular distribution was chosen based on the differential 
cross section distribution, as mentioned above. The 4- 
momentum of these events were then defined in the CM 
and boosted to the lab frame. With this procedure we 
determined the lab angular (0i a b) distribution. Then the 
neutrino-Earth interaction point (R) was chosen based on 
the interaction probability distribution. The separation 
at the detector is then simply 6i a bX R. 

The result is shown in Fig. [7] shows the NLSPs track 
separation distribution for various squark masses. The 
rather broad distribution is characterized by track sepa- 
rations of hundreds of meters. For 300 GeV squark mass, 
87% of the distribution is above 30 m separation, and 
52% is above 100 m separation. 

Another feature of the NLSP signal is that the two 
tracks should have very small relative angles. The angu- 
lar separation of the tracks in the signal is always well 
below a degree, the typical angular resolution of neutrino 
telescopes, so that the tracks will appear parallel to each 
other. Then, most NLSP events would consist of two par- 
allel but well separated tracks, and are therefore expected 
to be very distinctive and different from backgrounds, as 
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C. Backgrounds 
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FIG. 7: Top panel: Track separation distribution of In pair 
events. Bottom panel: The track separation distribution of 
the di-muon background. The relative normalization corre- 
sponds to the relative number of events for signal and back- 
ground. Note the different horizontal scales, as well as differ- 
ent binning between the two figures. 



The main potential background for the detection of 
NLSPs, muons produced by upgoing atmospheric neutri- 
nos, is eliminated by asking for two tracks in the events. 
With the additional requirement that the tracks be al- 
most parallel and considering the excellent time resolu- 
tion of the IceCube detector, the coincidence background 
(two coincident muons tracks) becomes negligible. 

However, there remains an important background, 
namely direct di-muon production. The main source of 
direct di-muons is the production of charm, which subse- 
quently decays semileptonically to a muon. That is, we 
consider the process 

vN — > (i~H c — > (i~ fi + H x v , 

where the charm hadron H c decays according to H c — > 
H x [i + v, and H x can be a strange or non-strange hadron. 
The cross section for producing a charm quark from a 
light d or s quark is 



d 2 a(vN -> cX) 
d£dy 



G\2ME V r2{q2h 



1 - 



2ME V £, 
|V^| 2 d(0} > 



{|v M | 2 *(0 

(4.28) 



where s(£) and d(£) are the strange and down quark par- 
ton distribution functions respectively. In eqn. (|4.28|) the 
slow scaling variable is defined as 



(4.29) 



and we also defined 



we show in Section [TVCj Given that the track separation 
ranges up to (300-400)m, and that the typical dimension 
of the detector is 1 Km, we expect that many events will 
be contained and therefore distinguishable from single- 
muon tracks. On the other hand, the typical separation 
appears to be large enough to distinguish the presence of 
two tracks in a given signal event. 



The energy with which the slepton pairs arrive at 
the detector has a distribution which peaks just above 
100 TeV, as it can be seen in Figure [5] On the other 
hand, the energy deposited in the detector by a typical 
stau track (~ 800 m which is IceCube geometrical effi- 
ciency for long tracks [l(3]) is approximately 160 GeV, 
which is small compared to its total energy. Low energy 
muons may mimic this energy deposition pattern, but 
they tend to have shorter paths in the detector. This 
fact can be used in order to reduce the background, as 
we discuss below. 



R{Q 2 ) 



M 2 W 



(4.30) 



The di-muon background results in a number of events 
larger than the signal, even for the case of a 300 GeV 
squark mass, as it can be seen in Table |U Also, the 
angular separation between the di-muon tracks is, as in 
the case of the NLSP signal, very small and well below 
the projected angular resolution of neutrino telescopes. 
Thus, di-muon tracks will also appear to be parallel. 



However, muons need to be close to the detector in 
order to range into it. Thus, the di-muon events have a 
considerably smaller separation between the two tracks 
than the NLSP events, which can range in from hun- 
dreds or even thousands of kilometers. This is shown 
in Figure [7l where it is clear that for track separations 
above 100 m there should not be any significant contri- 
bution from the di-muon background. We can estimate 
the statistical significance of the signal as S/y/~B, with 
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FIG. 8: Arrival energy distribution of the £r at the detector 
and for m- q = 300, 600 and 900 GeV.. Here, m- t = 150 GeV 
and ma — 250 GeV. Also shown is the arrival distribution 
for the di-muon background. The energy deposited in the 
detector by a stau traveling the average track length of 800 m 
is E^ cp — 160 GeV, approximately the same for all the masses 

considered here (See Figure [9}. 



S and B the number of signal and background events, 
respectively. We could then ask what is the cut in track 
separation that would yield a significance equal to 5. We 
find that cutting events with less than 106 m separation 
results in 3 detectable NLSPs and 0.25 di-muons, yielding 
a 5 er significance. 

It is possible to further reduce the di-muon background 
by making use of the energy deposition of the events. 
This is due to the fact that NLSPs go through the de- 
tector loosing very little energy, while muons tend to de- 
posit larger amounts of energy in the detector in such 
long paths. This can be seen by comparing the energy 
loss for a muon in Figure [3] with that of a slepton in 
Figure El where the largest source of energy loss for the 
NLSP is still three orders of magnitude smaller than that 
for the muon at the energies with which they arrive at 
the detector. For example, a muon with an arrival en- 
ergy of 3 x 10 3 GeV, transversing 800 m of ice (the av- 
erage through-detector length) will deposit ~ 450 GeV 
in its path. This is to be compared with the deposited 
energy for a slepton in Figure [9] Thus, excluding events 
with deposited energies above a few hundred GeV, we can 
eliminate the high energy tail of the di-muon background 
of Figure [5J This can be used, in addition to the track 



FIG. 9: Energy deposited in the detector by a right-handed 
stau traveling the average track length of 800 m. 



separation cut, to further reduce the background events. 
The precise extent to which this method can be used to 
reduce the background requires further study. 

Finally, there is an additional potential source of di- 
muons from v^j — > W + /i scattering [201 ] . where the 7 is 
emitted off the parton in the nucleon and the produced 
W decays to ijlv^. However, not only this is a smaller 
source of di-muons, but also its track separation distribu- 
tion is very similar to that of the di-muons coming from 
charm production and is equally eliminated by the same 
cut in track separation. 



V. CONCLUSIONS 

We have confirmed that the proposal of Ref . [l| , is vi- 
able, i.e. neutrino telescopes are potentially sensitive to 
the relatively long-lived charged NLSPs which are present 
in a wide variety of models of supersymmetry breaking. 
Compared to Ref. [l[ , we have now made a detailed study 
of the NLSP energy loss (Section 13 .4[) . This has con- 
firmed our estimates in Ref. [l[ and is in agreement with 
the calculations of Ref. Q . 

Regarding the analysis of the signal, we have shown 
that the separation between tracks will be such that we 
can expect that most events will be contained in a km 3 
such as IceCube, whereas they will be separated enough 
to be identified as two distinct tracks. Moreover, we have 
shown that the separation between tracks is a variable 
that will allow to identify the signal above the most im- 
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portant background, the direct production of pairs. 
This is shown in Figure [7] For instance, we have seem 
that a 5 a significance can be achieved with a cut of 106 m 
in the track separation. 

Also, given that the NLSP deposits energy in the de- 
tector like a rather low energy muon, but it goes right 
through it like a high energy muon could, it is possi- 
ble to further reduce the number of background events 
by making a cut in the amount of deposited energy, of 
the events going through the detector. The details of 
this method of reducing the background deserves further 
study. 

The region of the supersymmetry breaking parame- 
ter space that is available to neutrino telescopes is de- 
termined by the requirement that the NLSP lifetime be 
long enough to give a signal (V~F > 5 x 10 6 GeV). Thus 
the observation of NLSP events at neutrino telescopes 
will constitute a direct probe of the scale of supersym- 
metry breaking. This is to be compared with the po- 
tential observation of these NLSPs at the Large Hadron 
Collider (LHC), where for this range of \/F, the NLSP 
decays outside the detector and is seen through its ioniza- 
tion tracks. However, the observation at the LHC would 
not constrain the NLSP lifetime significantly. Thus, we 



see that neutrino telescopes are complementary to col- 
lider searches. For instance, the observation of NLSP 
events at the LHC, coupled to significantly fewer than 
expected events in neutrino telescopes would point to 
VF < 10 7 GeV. 

The event rates shown in Table U are already en- 
couraging for experimental facilities that are been built, 
such as IceCube. Future upgrades of IceCube [2l| will 
result in even better sensitivity. The water detectors 
ANTARES [12; NESTOR [H and NEMO [p will be 
also be able to look for these NLSPs. 

Finally, in the present work we focused on supersym- 
metry. However, other theories give rise to relatively long 
lived charg ed p articles which can be observed by neutrino 
telescopes [25[ . A large portion of this work could be ap- 
plied to these other scenarios with little modification. 
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